First Order Phase Transition and Phase Coexistence in a Spin-Glass Model 
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We study the mean-field static solution of the Blume-Emery-Griffiths-Capel model with quenched 
disorder, an Ising-spin lattice gas with quenched random magnetic interaction. The thermodynamics 
is worked out in the Full Replica Symmetry Breaking scheme. The model exhibits a high tempera- 
ture/low density paramagnetic phase. When the temperature is decreased or the density increased, 
the system undergoes a phase transition to a Full Replica Symmetry Breaking spin-glass phase. 
The nature of the transition can be either of the second order (like in the Sherrington-Kirkpatrick 
model) or, at temperature below a given critical value (tricritical point), of the first order in the 
Ehrenfest sense, with a discontinuous jump of the order parameter and a latent heat. In this last 
case coexistence of phases occurs. 

PACS numbers: 64.40.-i, 64.60.Cn, 75.10.Nr 



The spin-glass phase (SG) has played and is still play- 
ing a central role in the understanding of disordered and 
complex systems. The analysis of mean-field models re- 
vealed different possible scenarios for the SG phase and 
the transition to it. Most of the work, however, has 
been concentrated on just two of them. In order of 
appearance in literature the first scenario is described 
by a Full Replica Symmetry Breaking (FRSB) solution 
characterized by a continuous order parameter function 
, which continuously grows from zero by crossing the 
transition. The prototype model is the Sherrington- 
Kirkpatrick (SK) model ||, a fully connected Ising-spin 
with quenched random magnetic interactions. The sec- 
ond scenario, initially introduced by Derrida ||, provides 
a transition with a jump in the order parameter to a SG 
One Step Replica Symmetry Breaking (1RSB). No dis- 
continuity appear, however, in the thermodynamic func- 
tions. Actually, at the transition to the 1RSB SG phase 
the Edwards- Anderson order parameter can either grow 
continuously from zero or jump discontinuously to a finite 
value. The first case includes Potts-glasses with three of 
four states |Q, the spherical p-spin spin glass model in 
strong magnetic field || and some spherical p-spin spin 
glass model with a mixture of p = 2 and p > 3 inter- 
actions §, §. The latter case includes, instead, Potts- 
glasses with more than four states Q, quadrupolar glass 
models ft M, p-spin interaction spin-glass models with 
p > 2 ||, and the spherical p-spin spin glass model 

in weak magnetic field || . The models of this class, of- 
ten referred to as "discontinuous spin glasses" j[lj , have 
been widely investigated in the last years because of their 
relevance for the structural glass transition observed in 
fragile glasses Jl(| [l^] . 

In all cases discussed so far the transition is always 
continuous in the Ehrenfest sense. To our knowledge, 
the first case of a spin glass undergoing a genuine first 
order thermodynamic transition is the so called Ghatak- 
Sherrington model (GS) However, besides the full 



analysis of the RS solution, the study of the FRSB solu- 
tion for this model was done only close to the continuous 
transition (SK-likc) down to and including the tricritical 
point because of its complexity jl3[ [t4| . We also recall 
that an exactly solvable model, generalization of the Der- 
rida's REM ||, displaying a first order phase transition 
to a SG phase with latent head, was introduced by Mot- 
tishaw [^5[. However, at difference with the GS model, 
the SG phase is 1RSB. 

Recently a generalization of the GS model Jp| has been 
considered in connection with the structural glass tran- 
sition due to the conjectured existence Jl6| of a "discon- 
tinuous" , in the above mentioned sense, transition to a 
IRSB SG phase. This possibility has risen new inter- 
est in such model and its finite dimensions version has 
been numerically investigated in a search for evidence of 
a structural glass transition scenario |]l7f . 

To clarify this issue and its compatibility with previous 
results on the GS model, we have investigated the whole 
phase diagram, deep in the SG phase, for the mean field 
quenched disorder variant of the Blumc-Emery-Griffiths- 
Capel model jl8| (BEGC), introduced for the A transition 
in mixtures of He 3 -He 4 , which includes the GS model. 

In literature there exists two different versions of the 
model. The first one is the direct generalization of the 
original BEGC model and uses spin-1 variables <Ti = 
— 1, 0, 1 on each site i of a lattice [|l9[ ^0|, while the second 
formulation is a lattice gas (rii = 0, 1) of spin-1/2 vari- 
ables (Si = —1,1) 0, ^l). In both cases the spin vari- 
ables interact through quenched random couplings. The 
two formulations are equivalent, at least as far as static 
properties are concerned. Indeed, by imposing Ui = Si rii 
the two models can be transformed one into the other, 
apart from a rescaling of the chemical potential/crystal 
field 123] . In this paper we will use the second formulation 
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described by the Hamiltonian (DBEGC) pi 



n = - JijSiSjmrij - — mnj -nY^m (l) 



N 



*<3 



representing an Ising spin glass lattice gas coupled to a 
spin reservoir. The symmetric couplings Jy are quenched 
Gaussian random variables of zero mean and variance 
jfj = J?} = J 2 /N. Here and in the following the overline 
denotes average with respect to disorder. Limiting cases 
of the model are the SK model Q obtained for p,/ J — ► 00, 
the site frustrated percolation model |23| for K = —1 and 
J/fi — > 00, and the GS model for K = 0. 

To keep the level of the presentation as general as pos- 
sible, we shall avoid technical details and report only the 
main results for the phase diagrams for the three rele- 
vant cases of the GS model (K = 0), the frustrated Ising 
lattice gas (K = —1) and the case of attracting particle- 
particle interaction (K = 1). 

Applying the standard replica method, the FRSB solu- 
tion in the SG phase is described by the order parameter 
function IBJ 



q(x) = / dyP(x,y)m(x,y) 2 



and the density of occupied sites by 

cosh Jy 



f cc 
P = dyP(l,y) — 
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where &i = ((3 J) 2 [p - q(l)]/2 + /3(jU + Kp) and m(x, y) 
and P(x,y) p4] are solution of 

rn{x. y) = --m"(x, y) + A(x) m(x, y) m'(x, y) , (4) 



P(x,y) 



q(x) 



P"(x,y) + &(x)[P(x,y)m{x^ , (5) 



with boundary conditions 
m(l,y) = 



sinh(/3j/) 



e" 01 + cosh(/3y) ' 



P(0,y) = 



exp 



2?(0)J 



(G) 
(7) 



The functions m(x,y) and P(x,y) are, respectively, the 
local magnetization and local field probability distribu- 
tion at "time scale" x £ [0,1] j24j, while A(x) is the 
Sompolinsky's anomaly p5[ . As usual, the "dot" denotes 
partial derivative with respect to x while the "prime" the 
one with respect to y. 

All thermodynamic quantities can be written in terms 
of the above functions. For example, defining K = 
K + pj/2, the internal energy density u and the entropy 
density s read 



K 2 



H p- 



dx q(x) A{x) , (8) 



-pQi 



[p-s(i)r 



dyP(l,y) 



x {log [2 + 2 e 01 cosh(/3 Jy)] - (3Jym(l,y)} (9) 

We have solved the coupled equations (H)-© in the 
Parisi's gauge A = —(3Jxq(x) using the pseudo-spectral 
method introduced in |26| . Analyzing the stability of the 
RS solution one gets the critical lines 



1 - {pjpf = , 

1 - pk{\ - p) P = 



(10) 

(ii) 



above which the only solution is the paramagnetic (PM) 
solution q(x) = for x £ [0, 1], p = 1/[1 + e _Ql ]. This is 
stable for any value of K. In the T — p plane, these are, 
respectively, the straight line and the left branch of the 
spinodal line shown, in Figs. [TJ, ^ and|| (for K = 1, 0, — 1, 
respectively). The two lines meet at the tricritical point 



T c = Pc =- 3 / 2 + X+ ^ 2 -^ + 9 / 4 , (12) 



2K 



Ik: 



K + log 1 

Pc 



(13) 




FIG. l:T-p phase diagram of the DBEGC for K = 1. The 
dot marks the tricritical point p c = — 1, T c — 1/2, p c — 1/2. 
See text for discussion. 

By crossing the critical line ( |To| ) above the tricritical 
point (p > Pc, T > T c , fj, > fi c ) the system undergoes 
a continuous phase transition of the SK-type to a FRSB 
SG phase, with a non-trivial continuous order parameter 
function q{x) which smoothly grows from zero. 

Below the tricritical point the scenario is completely 
different with a transition from the PM phase to a FRSB 
SG phase with q{x) which discontinuously jumps from 
zero to a non-trivial (continuous) function. At the criti- 
cal temperature the entropy is discontinuous, see Fig. |^, 
and hence a latent heat is involved in the transforma- 
tion, implying that the transition is of the first order in 
the Ehrenfest sense. The transition line is determined 
by the free energy balance between the PM and the SG 
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FIG. 2: Entropy density as a function of temperature for 
K — 1. For p < fjb c = — 1 the entropy is discontinuous at the 
transition temperature. 
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FIG. 3: /i- p phase diagram of the DBEGC for K = 1. 
Three isothermal lines are plotted, two above and one below 
the tricritical temperature T c = 1/2. For T — 0.3 also the 
metastable branches are shown, both in the RS PM phase and 
in the FRSB SG phase. They reach the spinodal lines with 
zero derivative. In this plane of conjugated thermodynamic 
variables a Maxwell construction can be explicitly performed. 
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FIG. 4: T - p phase diagram of the DBEGC for K = 1. 



phase 1 15 , and is shown as a broken line in the phase 
diagrams. The line (|Tl|) where the PM solution becomes 
unstable, and the equivalent line from the SG side are 
the spinodal lines. This can be better appreciated in the 
fj, — p plane. From Fig. |^ we indeed see that the isother- 




FIG. 5: fj, - p phase diagram of the DBEGC for K = 0. 
The dot marks the tricritical point fi c — —0.731, T c = 1/3, 
p c = 1/3. The isothermal at T = 0.27 is plotted, together 
with its metastable parts (dotted line) both in the SG and in 
the PM phase. 
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FIG. 6: T — p and T — p. phase diagrams for K — 0. A line 
at constant p — —0.75 < p c is shown in the T — p plane. 
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FIG. 7: fi-p phase diagram of the DBEGC for K = -1, with 
isothermal at T = 0.2 < T c . The dot marks the tricritical 
point p c = -0.559, T c = 0.219, p c = 0.219. 

mal lines cross the instability lines with zero derivative 
and hence a diverging compressibility k — (1/ p )dp/dp,. 
It can be shown that the first order transition line can be 
determined in the p, — p phase diagram from the isother- 
mal and spinodal lines by using a Maxwell construction. 
In the region between the first order transition line and 
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FIG. 8: T — p and T — p phase diagrams for K = — 1. In the 
T — p plane the line at constant p = —0.57 < p c is plotted. 

the spinodal line the pure phase is metastable. Below 
the spinodal lines (in the T — p plane) no pure phase 
can exist and the system is in a mixture of PM and SG 
phase (phase coexistence). Finally, the phase diagram in 
the T — \i plane, for K = 1, is shown in Fig. |]. 

By varying K the scenario remains qualitatively un- 
changed. The only effect of a strong repulsive particle- 
particle interaction is to increase the phase diagram zone 
where the empty system (p = 0) is the only stable so- 
lution. In order to find further phases, e.g. an anti- 
quadrupolar phase Jig] , a generalization of the present 
analysis to a two component magnetic model p7|, in- 
cluding quenched disorder, has to be carried out p8fl. In 
Figs. we show the phase diagrams for K — 0, the GS 
model Jlq ], and K = — 1 the frustrated lattice gas p7J . 

In conclusion, we have discussed the complete phase 
diagram of the DBEGC model in the mean field limit, 
explicitly solving the FRSB equations in the whole SG 
phase with the pseudo-spectral method developed in Ref. 
[p6|i . Our results rule out the possibility of a 1RSB phase: 
the SG phase is always of FRSB type. The transition 
between the PM phase and the SG phase can be either 
of the SK-type or, below the tricritical temperature, a 
first order thermodynamic phase transition. In the latter 
case, like in the gas-liquid transition, a latent heat is 
involved in the transformation. Moreover, for a certain 
range of parameters (between the spinodal lines), no pure 
phase is achievable, not even as a metastable one, and the 
two phases coexist. 
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